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Abstract 

The giant planets of our solar system possess envelopes consisting mainly of hydro- 
gen and helium but are also significantly enriched in heavier elements relatively to 
our Sun. In order to better constrain how these heavy elements have been delivered, 
we quantify the amount accreted during the so-called "late heavy bombardment", 
at a time when planets were fully formed and planetesimals could not sink deep 
into the planets. On the basis of the "Nice model", we obtain accreted masses (in 
terrestrial units) equal to 0.15 ± 0.04 for Jupiter, and 0.08 ±0.01 M® for Saturn. 
For the two other giant planets, the results are found to depend mostly on whether 
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they switched position during the instability phase. For Uranus, the accreted mass 
is 0.051 ± 0.003 with an inversion and 0.030 ± 0.001 without an inversion. 
Neptune accretes 0.048 ± 0.015 in models in which it is initially closer to the 
Sun than Uranus, and 0.066 ± 0.006 otherwise. With well-mixed envelopes, this 
corresponds to an increase in the enrichment over the solar value of 0.033 ± 0.001 
and 0.074 ± 0.007 for Jupiter and Saturn, respectively. For the two other planets, 
we find the enrichments to be 2.1 ± 1.4 (w/ inversion) or 1.2 ± 0.7 (w/o inversion) 
for Uranus, and 2.0 ± 1.2 (w/ inversion) or 2.7 ± 1.6 (w/o inversion) for Neptune. 
This is clearly insufficient to explain the inferred enrichments of ~ 4 for Jupiter, 
~ 7 for Saturn and ~ 45 for Uranus and Neptune. 
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i 1 Introduction 



2 The four giant planets of our solar system have hydrogen and helium envelopes 

3 which are enriched in heavy elements with respect to the solar composition. 

4 In Jupiter, for which precise measurements from the Galileo probe are avail- 

5 able, C, N, S, Ar, Kr, Xe are all found to be enriched compared to the solar 



e value by factors 2 to 4 (jOwen et al. 



1999 



abundances based on the compilation by 



Wo ng et al. 



ratio is found to be 7.4 ±1.7 times solar (IFlasar et al. 



9 Neptune it is approximately 45 ± 20 times solar (iGuillot and Gautier . 



20041 ) (Assuming solar 



Loddersl (2003)). In Saturn, the C/H 



2005). In Uranus and 



20071) 
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(corresponding to about 30 times solar with the old solar abundances). In- 
terior models fitting the measured gravitational fields constrain enrichments 
to be between 1.5 and 8 for Jupiter and be tween 1.5 and 7 times the solar 



13 value for Saturn (jSaumon and Guillot 



2004 ) . For Uranus and Neptune, the 



envelopes are not massive enough (1 to 4 Earth masses) for interior models to 
provide global constraints on their compositions, 
rrriching giant planets in heavy elements is not straightforward. 



Guillot and Gladman 



20001 ) have shown that once the planets have their final masses, the ability of 



Jupiter to eject planetesimals severely limits the fraction that can be accreted 
by any planet in the system. The explanations put forward then generally 
imply an early enrichment mechanism: 



22 • 



Alibert et al.l (120051 ) show that migrating protoplanets can have access to 
a relatively large reservoir of planetesimals and accrete them in an early 
phase before they have reached their final masses and started their contrac- 
tion. This requires the elements to be mixed upward efficiently, which is 



26 energ etically possib 



27 core (IGuillot et al. 



e, an d may even lead to an erosion of Jupiter's central 



200J). 



28 • 



The forming giant planets may accrete a gas that has been enriched in heavy 
elements through the photoevaporation of the p rotoplanetary disk' s atmo- 



20061 ) 



30 sphere, mainly made of hydrogen and helium (IGuillot and Hueso 

31 This could explain the budget in noble gases seen in Jupiter's atmosphere 

32 but is not sufficient to explain the enrichment in elements such as C, N, 

33 O because small grains are preven ted from reach i ng th e planet due to the 



34 formation of a dust-free gap (eg. 



Paardekooper 



20071 ). The photoevapo- 



35 ration model requires that the giant planets form late in the evolution of 
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36 disk s, which app e ars consistent wi t 



see 



Ida and Linl ([2004]); 



Ida et al. 



i mo dern scenarios of planet formation 



( 120081 )). It also implies that a significant 



38 amount of solids are retained in the disk up to these lat e stages, as plausible 



from simulations of disk evolution (e.g. 



Garaudl f l2007h ) 



It has been recently suggested that the Solar System underwent a major 
change of structure during the phase calle d 'Late Heavy Bombardment' (LHB) 



( ITsiganis et al 



2005 



Gomes et al. 



20051 ). This phase, which occurred ~ 650 My 



after planet formation, was characterized by a spike in the cratering history 
of the terrestrial planets. 

The model that describes these structural changes, often called the 'Nice 
model' because it was developed in the city of Nice, reproduces most of the 
current orbital characteristics of both planets and small bodies. This model 
provides relatively tight constraints on both the location of the planets and on 
the position and mass of the planetesimal disk at the time of the disappear- 
ance of the proto-planetary nebula. Thus, it is interesting to study the amount 
of mass accreted by the planets at the time of the LHB, in the framework of 
this model. 

The article is organized as follows: we first describe the orbital evolution model 
at the base of the calculation. Physical radii of the giant planets at the time 
of the late heavy bombardment are also discussed. We then present results, 
both in terms of a global enrichment and in the unlikely case of an imperfect 
mixing of the giant planets envelopes. 
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5 8 2 The 'Nice' model of the LHB 



59 2.1 General description 



The Nice model postulates that the ratio of the orbital periods of Saturn and 
Jupiter was initially slightly less than 2, so that the planets were close to 
their mutual 1:2 mean motion resonance (MMR); Uranus and Neptune were 
supposedly orbiting the Sun a few AUs beyond the gas giants, and a massive 
planetesimal disk was extending from 15.5 AU, that is about 1.5 AU beyond 
the last planet, up to 30-35 AU. 

As a consequence of the interaction of the planets with the planetesimal disk, 
the giant planets suffered or bital migration, whi ch slowly increased their or- 



bital separation. As shown in 



Gomes et al. 



( 120051 ) N-body simulations, after a 



long quiescent phase (with a duration varying from 300 My to 1 Gy, depend- 
ing on the exact initial conditions), Jupiter and Saturn were forced to cross 
their mutual 1:2 MMR. This event excited their orbital eccentricities to values 
similar to those presently observed. 

The acquisition of eccentricity by both gas giants destabilized Uranus and 
Neptune. Their orbits became very eccentric, so that they penetrated deep 
into the planetesimal disk. Thus, the planetesimal disk was dispersed, and the 
interaction between planets and planetesimals finally parked all four planets 
on orbits with separations, eccentricities and inclinations similar to what we 
currently observe. 



79 This model has a long list o 



tecture of the giant planets ( iTsiganis et al. 



successes. It explain s the current orbital archi- 



20051 ) . It also explains the origin 



bi and the properties of the LHB. In the Nice model, the LHB is triggered by 
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the dispersion of the planetesimal disk; both the timing, the duration and the 
intensity of the LHB deduced from Lunar constraints are well reproduced by 



B4 the model ( iGomes et al. 



20051 ). 



Furthermore, the Nice model also explains the capture of planetesimals around 
the Lagrangian points of Jupiter, with a total mass and orbita l distr ibution 



20051 ). More 



consistent with the observed Jupiter Trojans ( IMorbidelfi et all 
recently, it has been shown to provide a framework for understanding the cap- 
ture and orbital distribution of the irregular s atelli tes of Saturn, Uranus and 



90 Neptune, but not of Jupiter ( iNesvorny et al. 



20071 ) . except in the case of an 



encounter between Jupiter and Neptune which is a rare but not impossible 
event. The main properties of the Kuiper belt (the relic of the primitive trans- 
planetary pl anetesimal disk) hav e als o been explained in the context of the 



94 Nice model (ILevison et al. 



(120081 ): see 



Morbidelli et al. 



( 120081 ) . for a review). 



2.2 The dynamical simulations 



96 In th is work, we use 5 of the numerical simulations performed in (I Gomes et al 



97 120051 ). The main simulation is the one illustrated in the figures of the Gomes 

98 et al. paper. The 1:2 MMR crossing between Jupiter and Saturn occurs after 

99 880 My, relatively close to the observed timing of the LHB (650 My). When 

100 the instability occurs, the disk of planetesimals still contained 24 of its initial 

101 35 Earth masses (M©). 

102 During the evolution that followed the resonance crossing, Uranus and Nep- 

103 tune switched position. Thus, according to this simulation, the planet that 

104 ended up at ~ 30 AU (Neptune) had to form closer to the Sun than the 

105 planet that reached a final orbit at ~ 20 AU (Uranus). 
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106 However, because the planets evolutions are chaotic during the instability 

107 phase, different outcomes can be possible. Thus Gomes et al. performed 4 
los additional simulations with initial conditions taken from the state of the sys- 
109 tern in the main simulation just before the 1:2 resonance crossing, with slight 
no changes in the planets' velocities. Two of these 'cloned' simulations again 
in showed a switch in positions between Uranus and Neptune, but the two oth- 
n2 ers did not. That is, in these two cases the planet that terminated its evolution 
in at 30 AU also started the furthest from the Sun. 

H4 Here we use these 5 simulations (the main one and its 4 'clones') to evaluate 

us the amount of solid material accreted by the planets and how it could vary 

lie depending on the specific evolutions of the ice giants. Notice that, whereas the 

in main simulation spans 1.2 Gy (and therefore continues for 320 My after the 

us 1:2 MMR crossing), the cloned simulations cover only a time-span of approx- 

n9 imately 20 My, and were stopped when the planets reached well separated, 

120 relatively stable orbits. 

121 2.3 Probability of impact and accretion of planetesimals 

122 From these dynamical simulations, several steps are necessary to estimate the 

123 amount of mass accreted by the planets. 

124 For each simulation at each output time (every 1 My) we have the orbital 

125 elements of the planets and of all the planetesimals in the system. We are 

126 aware that this time interval may be a bit too long to resolve the evolution 

127 of the system during the transient phases that immediately follow the onset 

128 of the planetary instability. On the other hand, when the instability occurs, 

129 most of the disk is still located beyond the orbits of the planet, so that the 
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130 bombardment rate is not very high. Thus, we believe that this coarse time 

131 sampling is enough for our purposes. 

132 First we look for planetesimals that are in a mean motion resonance with a 

133 planet. The resonances taken into account are the 1:1, 1:2, 2:3, 2:1, 3:2. When 

134 computing the collision probability with a planet, the objects in resonance 

135 with that planet will not be taken into account (but they will be considered 

136 for the collision probability with the other planets). The rationale for this 

137 is that the resonant objects, even if planet-crosser, cannot collide with the 

138 planet, because they are phase-protected by the resonant configuration, as in 

139 the case of Pluto. 



140 The width of a resonance is proportionnal to y ^' s anct * a where M p i anet is 

hi the mass of the considered planet and a is the semi-major axis of the precise 

142 resonant orbit. Hence we take an approximative relative width of — = 1% on 

143 the semi-major axis to define the area where a planet and a given planetesimal 

144 are considered to be in resonance. Then we select all the non-resonant particles 

145 that cross the orbit of a planet. The intrinsic collision probability Pj of each 



146 of these particle s with the planet is computed using t he method deta i 



Farinella et al. 



ed in 



(1992) 



147 IWetherilll ( 119671 ). implemented in a code developed by 

148 and kindly provided to us. Once Pj is known for each particle (i = 1, . . . , N) 

149 the mass accreted by the planet during the time-step At is simply : 



N 



151 where Mj is the mass of the planetesimal (with Mj = 0.00349 Earth mass), 

152 and f gr av is the gravitational focusing factor. The latter is equal to : 

153 fgrav = 1 + TFTf (2) 

^rel 
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154 where V ie \ is the relative velocity between the planet and the planetesimal, and 

155 Viib is the escape velocity from the planet. Finally, the total mass accreted by 
we a planet during the full dynamical evolution is simply the sum of M acc over 

157 all time-steps taken in the simulation. 

158 3 Results: mass accreted by each giant planets 

159 Fig.l shows the cumulated mass captured by Jupiter, Saturn, Uranus and 
wo Neptune as a function of time in the case of the main simulation. The abrupt 
lei increase at 882 Ma is due to the triggering of the LHB when Saturn crosses 

162 the 1:2 resonance with Jupiter. It is interesting to notice that this short phase 

163 accounts for about two third of the mass acquired by the planets during their 

164 full evolution. 

165 Qualitatively, the more massive is the planet, the larger is the massed accreted 
lee from the planetesimal disk. This is because larger planets have larger gravita- 

167 tional cross sections. 

168 Uranus and Neptune have comparable masses, and therefore which planets 

169 accretes more mass depends on their orbital histories. In the model shown 
no in Fig.l, Uranus first accretes planetesimals at a larger rate than Neptune 

171 because Uranus is initially the furthest planet in the system and the closest 

172 to the planetesimal disk. However, when the two planets exchange position 

173 and Neptune is scattered into the planetesimal disk, it accretes many more 

174 planetesimals and eventually exceeds Uranus in terms of total accreted mass. 

175 As previously described, the evolution of the system is very chaotic and Fig.l 

176 only represents one of the possible outcomes. In order to assess the variability 

177 of the solutions, we also present in Fig. 2 the evolution of the 4 "cloned" sim- 
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178 ulations focused around the critical period of the 1:2 MMR crossing. We note 

179 them 'simu a\ 'simu b\ 'simu c', and 'simu d\ These simulations were started 

180 at 868 Myr, just before the 1:2 MMR crossing 880 Myr), and stopped at 

181 893, 897, 875 and 899 Myr respectively, as soon as the planets reached well 

182 separated and relatively stable orbits. 

183 Fig. 2 shows that the variability of the accreted masses during that period 

184 amounts to up to a factor 2 for all planets except Neptune, for which the 

185 variability is a factor 4. The added uncertainty on the results due to the 1 

186 Myr timestep appears small in comparison, as shown by the regularity of the 

187 curves. 

las In order to obtain the evolution of the mass accreted by 4 giant planets during 

189 the entire 1.2 Gyr period and to assess the effect of the position switch be- 

190 tween Uranus and Neptune on the final accreted mass, we proceed as follows: 

191 we simply assume that the planets accreted the same amount of mass as in 

192 the main simulation over the first 880 Myr and over the time ranging from the 

193 end of each cloned simulation up to 1,200 Myr. In the cases in which Uranus 

194 and Neptune do not switch positions, we consider that Uranus accreted before 

195 the LHB the same mass accreted by "Neptune" (the 4th planet) in the main 
we simulation, and inversely for Neptune. The results are shown in Fig. 3. We 

197 can see that the results for simulations are very similar to the result for the 

198 main simulation. However in the simulation a, Neptune eventually accretes 

199 less mass than Uranus. Conversely, the results for the simulations b and d are 

200 qualitatively different. Neptune is initially the closest planet to the disk and 

201 hence accretes much more planetesimals than Uranus also before the LHB. 

202 This remains the case during/after the LHB, since Neptune is scattered into 

203 the disk and acquires even more planetesimals compared to Uranus. 

204 Table 1 summarizes the total masses accreted by the planets, and compares 
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205 them to the masses of heavy elements in their hydrogen-helium envelopes es- 

206 tim ated from inte rior models fitting the giant planets gravitational moments 



207 sec 



Guillot 



20051 ). As before, for Uranus and Neptune, we separate the cases 

208 in which these planets exchange their positions from the cases in which they 

209 do not. 

210 In the first case, Uranus accretes an amount of planetesimals comparable to 
2n Neptune's, whereas when the order of the ice giants is not switched, Neptune 

212 accretes twice more planetesimals than Uranus. In all cases, the masses ac- 

213 creted are significantly smaller than the masses of the envelopes. 

214 For Jupiter and Saturn, the mass accreted is much lower 10~ 3 times 

215 smaller), whereas this ratio can increase to ~ 7 x 10~ 2 for Uranus and Nep- 

216 tune. Therefore in the framework of the 'Nice' model, the LHB has a stronger 

217 impact in terms of heavy elements supply relatively to the envelope mass, in 

218 the case of the two latter planets, Uranus and Neptune. 



219 4 Calculation of the envelope enrichments 



220 4-1 Fully mixed case 



221 Once we calculated the mass that each planet accreted during this period, 

222 it is straightforward to infer the corresponding change in composition. We 

223 thus calculate the increase of the atmosphere's enrichment AS, defined as the 

224 amount of heavy elements for a given mass of atmosphere compared to that 

225 same value in the Sun. More specifically, the global enrichment increase of 

226 a giant planet envelope of mass M enve i opo accreting a mass of planetesimals 
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227 -^^accreted 

(assuming that planetesimals don't reach the core) is: 



AS] 



M, 



accreted 



LHB 



M. 



envelope 



x Zp 



(3) 



229 where Z & is the mass fraction of heavy elements in the Sun. Following 



Grevesse et al 



230 ( 120051 ). we use in mean Z G = 0.015. This global enrichment is also the enrich- 

231 ment of the atmosphere, provided the envelope is well-mixed, a reasonnable 

232 ass umption given the fact that these planets should be mostly convective (see 



233 e.g. 



Guillot 



20051 ). These values of enrichment are calculated by taking the 

234 mean of the accreted masses of the table 1, and the uncertainty on the enve- 

235 lope mass is taken into account. Table 2 shows that this yields relatively small 

236 enrichments: the contribution of this late veneer of planetesimals accounts for 

237 only about 1% of the total enrichments of Jupiter and Saturn, and up to 10% 



238 in the case of Uranus and Neptune, owing to their smaller envelopes. 



4-2 Incomplete mixing case 



240 Mixing in the envelopes of giant planets is expected to be fast compared to the 

241 evolution timescales, and rather com plete because these planets are expected 



242 to be fully convective (j Guillot 



20051 ) . We want to test the possibility, however 



243 unlikely, that mixing was not complete, and that the observed atmospheric 

244 enrichments were indeed caused by these late impacts of planetesimals. 

245 The values of enrichment, in the hypothesis of an incomplete mixing of the 

246 envelope, depends on two elements : the extent of mixing of heavy elements in 

247 the envelope, but also the penetration depth of planetesimals in the envelope 

248 as a function of their size distribution at the time of the LHB. 

249 First let us evaluate to what extent the mixing should occur in order to retrieve 
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250 the observed enrichments (£c/h in Table 2). Following Eq.3, we evaluate the 

251 mass of envelope over which planetesimals should penetrate and be mixed to 

252 explain the observations : 

M-M-accrctcd , , 

mixed = 7 (,4J 

254 Using hydrostatic balance, and assuming a constant adiabatic gradient of 0.3, a 

255 pressure of 1 bar at the top of the atmosphere and a gravitational acceleration 

256 constant and equal to its value at the top of the atmosphere, we calculate the 

257 corresponding penetration depth, h mixc d, and pressure, P mixed, at this depth . 

258 We now have to consider in addition the penetration depth of planetesimals 

259 in the envelope as a function of their size distribution at the time of the LHB. 

260 The reason is that, even if we could define the extents of penetration and 

261 mixing giving the observed enrichments, an important fraction of planetesi- 

262 mals penetrating more deeply in the envelope would anyway cause a heavy 

263 elements supply over a larger extent, implying atmospheric enrichments lower 

264 than those observed. 

265 We thus have to determine the mass of the envelope shell enriched by a plan- 

266 etesimal of a given size. The main assumption here is to consider that during 

267 its entry into the atmosphere, the planetesimal desintegrates and mixes com- 

268 pletely with the atmosphere after crossing a mass of gas column equal to its 

269 own mass. Thanks to a parallel plane approximation of the atmosphere, the 

270 mass of the atmospheric shell thus enriched can be inferred from the ratio 

271 between the planetary area a p i anet and the planetesimal section a p \ multiplied 

272 by the mass of the considered planetesimal M p i : 



-^enriched shell — M p \ plaDet (5) 

a p i 
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274 We now define s m ; xe d, the critical planetesimal radius for which M enric h e d shell = 

275 M mixe d. Following Equations 5 and 6 and considering ice spherical planetesi- 

276 mals with a mass density noted p : 



■5 mixed 



3 X -Mmixed 
4 X Gtpianet X p 



(6) 



278 All the planetesimals larger than s m i xc d will penetrate more deeply than the 

279 extent of mixing and will enrich a larger part of the envelope, yielding the 

280 enrichments lower than observed. 

281 We now evaluate the mass fraction of planetesimals with sizes larger than 

282 s mixc d. For that we use a bi- modal size distribution inspir ed from the obser- 



283 vat ions of Trans- Neptunian Bodi e s ([Bernstein et al. 



284 fully by 



2004J), and used success- 



Charnoz and Morbidellil ( 120071 ) to explain the number of comets in 



285 the scattered disk and the Oort cloud : 



dN_ f 

i J small ^ $ 

as 

~ds~ = fbl9 X S ~ 4 ' 5 



-3.5 



S < S 
S > S 



(7) 



288 with 50km < sq < 100km the turnover radius, and f sma n and fbi g the nor- 

289 malisation factors whic h depends on the v alue of so and the total mass of 



290 the planetesimals disk (IGomes et al. 



20051 ). For each planet, the mass frac- 



291 tion is calculated with so=50 and 100 km in order to have a good range of 

292 values around the estimate d one which is approximately 70 km according to 



293 iFuentes and Holmanl (120081 ). 

294 Table 3 summarizes the results obtained in the case of an incomplete mixing. 

295 Compared to the whole envelope mass, the masses of layer enriched by this 

296 incomplete mixing are of the order of 1% for Jupiter and Saturn, and between 
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297 5 and 10% for Uranus and Neptune. As previoulsy mentioned, these values and 

298 those of the related quantities are a priori unrealistic because of the globally 

299 convective structure of the giant planets. Moreover, according to the results of 

300 mass fraction in Table 3, we see that the large planetesimals with a size bigger 

301 than s m ixcd are comparable and even predominant in terms of mass compared 

302 to the small ones, especially for Uranus and Neptune. 

303 Therefore even if we assume an incomplete mixing giving the observed en- 

304 richments, the important supply of heavy elements by the large planetesimals 

305 at layers deeper than h mixcd will imply anyway lower enrichments than those 

306 observed. 

307 In summary, it appears that the observed enrichments cannot be explained in 

308 the context of the Late Heavy Bombardment even by using the hypothesis of 

309 an incomplete mixing. 

310 5 Conclusion 

3n In this work, we evaluated the extent to which the late heavy bombardment 

312 could explain the observed enrichments of giant planets. 

313 We calculated the mass accreted by each planet during this period thanks to 

314 several dynamical simulations of the LHB within the so-called "Nice" model. 

315 The accreted masses were found to be much smaller than those of the en- 

316 velopes of each giant planet. In the realistic hypothesis of a global mixing in 

317 these envelopes, we found the enrichments over the solar value to be approxi- 

318 mately two orders of magnitude smaller than the observations for Jupiter and 

319 Saturn and one order of magnitude smaller than the observations for Uranus 

320 and Neptune. 
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321 We then tested the possibility of an incomplete mixing in the giant planets 

322 envelopes to account for the observed enrichments. With a size distribution of 

323 planetesimals inferred from observations of trans-neptunian bodies, we found 

324 that the enrichments were always at least a factor of 2 lower than observed. 

325 Given the efficient convection expected in the deep atmosphere, we expect 

326 however the mixing to be complete. 

327 Therefore we conclude that the enriched atmospheres of the giant planets do 

328 not result from the Nice model of the LHB and p robab ly from any model de- 



329 scribing the LHB. In fact 



Guillot and Gladmanl ( 120001 ) 's calculations showed 



330 that the mass needed to explain Jupiter's and Saturn's enrichments would 

331 be certainly much too large, in any late heavy bombardment model. Earlier 

332 events should thus be invoked in the explanation of the enriched atmopspheres 

333 of giant planets. On the other hand the enrichment process during the LHB 

334 may not be completely negligeable when considering fi ne measurements of the 



335 compositions of giant planets (eg. 



Marty et al 



20081 1 . When present it may 



336 also have a role in enriching the envelopes of close-in extrasolar giant planets 

337 because of their radiative structure. 



340 We acknowledge support from the Programme National de Planetologie. We 

341 thank one of the referees, Brett Gladman, for comments that improved the 

342 article. 



343 References 

344 Alibert, Y., Mordasini, C, Benz, W., Winisdoerffer, C, Apr. 2005. Models of 

16 



345 giant planet formation with migration and disc evolution. A&A 434, 343- 

346 353. 

347 Bernstein, G. M., Trilling, D. E., Allen, R. L., Brown, M. E., Holman, M., 

348 Malhotra, R., Sep. 2004. The Size Distribution of Trans- Neptunian Bodies. 
AJ 128, 1364-1390. 

350 Charnoz, S., Morbidelli, A., Jun. 2007. Coupling dynamical and collisional 

351 evolution of small bodies. Icarus 188, 468-480. 

352 Farinella, P., Davis, D. R., Paolicchi, P., Cellino, A., Zappala, V., Jan. 1992. 

353 Asteroid collisional evolution - an integrated model for the evolution of 

354 asteroid rotation rates. A&A 253, 604-614. 

355 Flasar, F. M., Achterberg, R. K., Conrath, B. J., Pearl, J. C, Bjoraker, G. L., 

356 Jennings, D. E., Romani, P. N., Simon-Miller, A. A., Kunde, V. G., Nixon, 

357 C. A., Bezard, B., Orton, G. S., Spilker, L. J., Spencer, J. R., Irwin, P. G. J., 

358 Teanby, N. A., Owen, T. C, Brasunas, J., Segura, M. E., Carlson, R. C, 

359 Mamoutkine, A., Gierasch, P. J., Schinder, P. J., Showalter, M. R., Fer- 

360 rari, C, Barucci, A., Courtin, R., Coustenis, A., Fouchet, T., Gautier, D., 

361 Lellouch, E., Marten, A., Prange, R., Strobel, D. F., Calcutt, S. B., Read, 

362 P. L., Taylor, F. W., Bowles, N., Samuelson, R. E., Abbas, M. M., Raulin, 

363 F., Ade, P., Edgington, S., Pilorz, S., Wallis, B., Wishnow, E. H., Feb. 2005. 

364 Temperatures, Winds, and Composition in the Saturnian System. Science 
307, 1247-1251. 

366 Fuentes, C. I., Holman, M. J., Jul. 2008. a SUBARU Archival Search for Faint 

367 Trans- Neptunian Objects. A J 136, 83-97. 

sea Garaud, P., Dec. 2007. Growth and Migration of Solids in Evolving Protostel- 

369 lar Disks. I. Methods and Analytical Tests. APJ 671, 2091-2114. 

370 Gomes, R., Levison, H. F., Tsiganis, K., Morbidelli, A., May 2005. Origin of 

371 the cataclysmic Late Heavy Bombardment period of the terrestrial planets. 



17 



372 Nature 435, 466-469. 

373 Grevesse, N., Asplund, M., Sauval, A. J., 2005. The New Solar Chemical Com- 

374 position. In: Alecian, G., Richard, O., Vauclair, S. (Eds.), EAS Publications 

375 Series. Vol. 17 of Engineering and Science, pp. 21-30. 

376 Guillot, T., Jan. 2005. THE INTERIORS OF GIANT PLANETS: Models and 

377 Outstanding Questions. Annual Review of Earth and Planetary Sciences 33, 

378 493-530. 

379 Guillot, T., Gautier, D., 2007. Giant planets. Treatise on Geophysics, Vol. 10, 

380 Planets and Moons, Elsevier, pp. 439-464. 

ssi Guillot, T., Gladman, B., 2000. Late Planetesimal Delivery and the Composi- 

382 tion of Giant Planets (Invited Review). In: Garzon, G., Eiroa, C, de Winter, 

383 D., Mahoney, T. J. (Eds.), Disks, Planetesimals, and Planets. Vol. 219 of 

384 Astronomical Society of the Pacific Conference Series, pp. 475 — h 

385 Guillot, T., Hueso, R., Mar. 2006. The composition of Jupiter: sign of a (rela- 
sse tively) late formation in a chemically evolved protosolar disc. MNRAS 367, 
387 L47-L51. 

sss Guillot, T., Stevenson, D. J., Hubbard, W. B., Saumon, D., 2004. The interior 

389 of Jupiter. Jupiter. The Planet, Satellites and Magnetosphere, pp. 35-57. 

390 Ida, S., Guillot, T., Morbidelli, A., Oct. 2008. Accretion and Destruction of 

391 Planetesimals in Turbulent Disks. APJ 686, 1292-1301. 

392 Ida, S., Lin, D. N. C, Mar. 2004. Toward a Deterministic Model of Planetary 

393 Formation. I. A Desert in the Mass and Semimajor Axis Distributions of 

394 Extrasolar Planets. APJ 604, 388-413. 

395 Levison, H. F., Morbidelli, A., Vanlaerhoven, C, Gomes, R., Tsiganis, K., 

396 Jul. 2008. Origin of the structure of the Kuiper belt during a dynamical 

397 instability in the orbits of Uranus and Neptune. Icarus 196, 258-273. 

398 Lodders, K., Jul. 2003. Solar System Abundances and Condensation Temper- 



18 



399 atures of the Elements. ApJ 591, 1220-1247. 

wo Marty, B., Guillot, T., Coustenis, A., Achilleos, N., Alibert, Y., Asmar, S., 

401 Atkinson, D., Atreya, S., Babasides, G., Baines, K., Balint, T., Banfield, D., 

402 Barber, S., Bezard, B., Bjoraker, G. L., Blanc, M., Bolton, S., Chanover, N., 

403 Charnoz, S., Chassefiere, E., Colwell, J. E., Deangelis, E., Dougherty, M., 

404 Drossart, P., Flasar, F. M., Fouchet, T., Frampton, R., Franchi, I., Gau- 

405 tier, D., Gurvits, L., Hueso, R., Kazeminejad, B., Krimigis, T., Jambon, 

406 A., Jones, G., Langevin, Y., Leese, M., Lellouch, E., Lunine, J., Milillo, 

407 A., Mahaffy, P., Mauk, B., Morse, A., Moreira, M., Moussas, X., Murray, 

408 C, Mueller- Wodarg, I., Owen, T. C, Pogrebenko, S., Prange, R., Read, 

409 P., Sanchez-Lavega, A., Sarda, P., Stam, D., Tinetti, G., Zarka, P., Zar- 

410 necki, J., Schmidt, J., Salo, H., Aug. 2008. Kronos: exploring the depths of 

411 Saturn with probes and remote sensing through an international mission. 

412 Experimental Astronomy, 34 — h 

413 Morbidelli, A., Levison, H. F., Gomes, R., 2008. The Dynamical Structure 

414 of the Kuiper Belt and Its Primordial Origin. The Solar System Beyond 
«5 Neptune, pp. 275-292. 

416 Morbidelli, A., Levison, H. F., Tsiganis, K., Gomes, R., May 2005. Chaotic 

417 capture of Jupiter's Trojan asteroids in the early Solar System. Nature 435, 

418 462-465. 

419 Nesvorny, D., Vokrouhlicky, D., Morbidelli, A., May 2007. Capture of Irregular 

420 Satellites during Planetary Encounters. AJ 133, 1962-1976. 

421 Owen, T., Mahaffy, P., Niemann, H. B., Atreya, S., Donahue, T., Bar-Nun, 

422 A., de Pater, I., Nov. 1999. A low-temperature origin for the planetesimals 

423 that formed Jupiter. Nature 402, 269-270. 

424 Paardekooper, S.-J., Jan. 2007. Dust accretion onto high-mass planets. A&A 

462, 355-369. 



19 



426 Saumon, D., Guillot, T., Jul. 2004. Shock Compression of Deuterium and the 

427 Interiors of Jupiter and Saturn. ApJ 609, 1170-1180. 

428 Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H. F., May 2005. Origin of 

429 the orbital architecture of the giant planets of the Solar System. Nature 

430 435, 459-461. 

431 Wetherill, G. W., 1967. Collisions in the asteroid belt. JGR 72, 2429-2444. 

432 Wong, M. H., Mahaffy, P. R., Atreya, S. K., Niemann, H. B., Owen, T. C, Sep. 

433 2004. Updated Galileo probe mass spectrometer measurements of carbon, 

434 oxygen, nitrogen, and sulfur on Jupiter. Icarus 171, 153-170. 

435 Caption figure 1 : Mass accreted (in Earth mass units) by Jupiter (plain), 

436 Saturn (dashed), Uranus (dash-dotted) and Neptune (dotted) respectively, as 

437 a function of time (in years). The simulation corresponds to the main simula- 

438 tion described in the text, in which Uranus and Neptune switch their relative 

439 positions. 

440 caption figure 2 : Additional mass accreted (in Earth mass units) during the 

441 time range of the 4 'cloned' simulations. Figures a and c (left) correspond 

442 to the cases in which Uranus and Neptune exchange position at the time of 

443 the LHB. Figures b and d (right) show the result of simulations in which the 

444 four planets preserve their initial order. These 'cloned' simulations start 868 

445 Myr after the beginning of the planets migration and are stopped once giant 

446 planets acquired well separated and stable orbits. 

447 Caption figure 3 : Mass accreted(in Earth mass units) for the 4 'cloned' sim- 

448 ulations during the whole time scale of the 'Nice' model. In each of these four 

449 panels, the period before 868 Myr and after 875-899 Myr (depending on the 

450 simulation) is assumed to be identical to the main simulation. Figures a and c 
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451 
452 
453 



(left) correspond to the cases in which Uranus and Neptune exchange position 
at the time of the LHB. Figures b and d (right) show the result of simulations 
in which the four planets preserve their initial order. 



454 Caption table 1 : Planetesimal masses accreted by the giant planets after the 

455 disappearance of the protosolar gaseous disk. 



456 Caption table 2 : Enrichment increase (see equation 3) calculated from the 
ation s of the model. The values of the observed ec/u are derived 



457 different simu 

458 from 



Guillot 



(120051 ) . 



459 Caption table 3 : Mass of enriched layer, extent of mixing h mixed , pressure 

460 level at the bottom of the mixing area P mixed, and planetesimal radii s m i xe d, 

461 which would match the observed enrichments. The last column is the mass 

462 percentage corresponding to the planetesimals whose the radius is larger than 

463 s m i X ed, the range being due to the two limiting values of so used. 
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Fig. 2. 
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Table 1 



Giant planet 


Envelope mass [M®] 


Accreted mass [M®] 


Jupiter 


300 - 318 


0.11 - 0.20 


Saturn 


70-85 


0.06-0.10 


Uranus (w/ inversion) 
(w/o inversion) 


1-4 


0.048 - 0.055 
0.029 - 0.031 


Neptune (w/ inversion) 
(w/o inversion) 


1-4 


0.033 - 0.064 
0.060 - 0.072 
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Fig. 3. 




Time (billion years) Time (billion years) 



Table 2 



Giant planet 


Ae L HB 


Observed £c/h 


Jupiter 


0.032-0.034 


4.1 ± 1 


Saturn 


0.062-0.076 


7.4 ± 1.7 


Uranus (w/ inversion) 
(w/o inversion) 


0.8-3.4 
0.5-2 


45 ± 20 


Neptune (w/ inversion) 
(w/o inversion) 


0.8-3.2 
1.1-4.4 


45 ± 20 
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Table 3 



Giant planet 


Mmixcd 


hmixcd 


P mixed 


Smixed 


% M (s > s mixcd ) 




(M ffi ) 


(km) 


(bar) 


(km) 




Jupiter 


3.41 


2000 


48300 


248 


23%-32% 


Saturn 


0.82 


2200 


7200 


86 


39%-54% 


Uranus (w/ inversion) 


0.11 


1200 


4200 


63 


44%-60% 


(w/o inversion) 


0.07 


1000 


2500 


37 


57%-69% 


Neptune (w/ inversion) 


0.11 


1000 


5100 


61 


45%-61% 


(w/o inversion) 


0.15 


1400 


13100 


84 


39%-54% 
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